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PROGRAM TITLE: RECOMP II BESSEL FUNCTION OP THE SECOND KIND, 

ORDER ONE SUBROUTINE (FLOATING POINT, RELOCATABLE) 

.1. PURPOSE 

1.1 To compute the Bessel function of the second kind, order one of a 
floating point argument, found in the A and R registers. 

2. METHOD 

2.1 For < X < 8, 

N 
Nl (x)« f (Gamm*Log e f )■ J^x) - |j - | £ A n > 

n=0 

where 

-12 
and N is the smallest integer for which A. T < 10- . . 

2o2 For 8 < x <oo 

1/2 



V X >=* [if] Un(X-f) 



N 



Z (=1)n [Vi^^- f > * h* Bin < x - f >] 



n«i 

where 

l 2 . 3 s . 5 2 ... (2 i -3) 2 (2 i -l)(2 i +l) 



B i 



i: (8X) 1 



and where N is the first integer for which 
whichever occurs first. 
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RESTRICTIONS 

3.1 Range of Argument X > 

3.2 Contents of L and V loops are destroyed, 

3.3 The following subroutines should be in memory: sin-cos(AN-012.l) 
log e (AN-037), J^AN-056). 

3. 1 * The following modifications should be made in AN-056 so that it will 
use AN-012.1 



L(0rigin plus 0011) 



k. USE 

k.1 Calling Sequence 



TRA 0000 
N0P 



N0P 
N#P 



a 


FCA 
TRA 


ARG 
N l 


m-1 


PZE 

PZE 


L{AN-CXL2.1) 
L(AN-056) 


Gtf-2 


PZE 
ERROR 


L(ATT-037) 
RETURN 


G#3 


NORMAL 


RETURN 



h.2 This routine occupies words 0000 thru OI67 or 120 full words. 

h.3 Minimum Accuracy: 8 decimal places. If speed is desired and several 

-P 

fewer accurate places are needed, it is suggested that 10 be entered 

into location 0i62~3^ where P is one more than the number of accurate 
decimal places required. 

k.k After a noxra&l exit the value of the function is found in the A and R 
registers. Error return is caused by negative argument. 
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5. CODING JNFOBtmHOU 
5.1 Constants 

5.1.1 Floating Point 



1.0 


in 


0002 


2.0 


in 


00<& 


8.0 


in 


OC0.O 


2/n « O.6366I 9772 


in 


0050 


Gamma « 0. 57721 5661*9 


in 


OO56 


3*A =* 


in 


00.61* 


3.0 


in 


0166 


10- 12 


in 


0162 



5.1.2 Fixed Point 

1 at B39 in 0003 
1 at B38 in 0005 

5.2 Erasable Locations: 0006-7, 005^-5, QI36-7, OI56-7. 

5.3 Master Hape includes: 

1) Basic routine, 0000-0167 

2) Relocation natrix 3777-^167 

3) AN- 00^ Relocation routine, 772*1-7757 • 
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Set Up Exit 
Addresses 



X : 



M 



JL 



X : 8 



>8 



y» 




a a 1 

SUM '« 



JL 



SUM - 



1/2 



±. 



C . » B rt _ cos (X - -^ ) 
n 2n-l * ■ TT ' 



♦ B^ sln< 



&) 

i~-; 



«r 




n=a+l 




^9 



»j SUM = jj(<Jamma + Los e fjJjOO-J 

3: 



n a 



Compute A 



SUM a SUM + A 




n=n+l 



A n :10 



,-12 



>10 



-IS 



<10 



rl2 




